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I. INTRODUCTION 

Charge exchange has been demonstrated to play an important role as 
a pumping mechanism for population inversions in a variety of near-thermal 
metal vapor lasers' 1 . The general process occurring is 

R + + M ~ — >■ R + M + * + A f (1) 

where R is a rare gas atom and M is a metal atom. The quantity Af is 
the energy defect of the state f for the metal ion (a^ > 0 corresponds 
to an increase in kinetic energy of the infinitely separated quasi - 
molecule). The rate for such a process is given by 

Rate (f) va f (2) 

where v is the relative velocity of the colliding pair and Of is the 

cross section for production of state f. One might expect that the cross 

section is largest for those states f with smallest ener».> / vfect Af, 

i.e. , for near resonant collisions. 

o 

Chubb and Rose have suggested that the rate of charge exchange will 
be enhanced if the relative collision velocity v is increased over thermal 
velocities by using a MPD arc accelerator as a high density source of 
rare gas ions. An important question then is whether the efficiency of 
the lasing system can be enhanced by assuring that the charge exchange 
takes place selectively to a few states f (ideally to a single state). 

The notion that near resonance can provide such a selection mechanism 
must be investigated more thoroughly. 
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As pointed out by Melius^, charge transfer occurs most efficiently 
in the region of atom-ion separation where the potential energies of the 
quasi -molecular states cross (or nearly cross) one another. Therefore, 
even though may be small there is no assurance that the potential 
curves cross in a region of strong coupling. Further there may be several 
final states f for v/hich curve crossings occur in a region of strong 
coupling. Thus a thorough understanding of the charge exchange process 
requires a more detailed analysis involving the relevant potential energy 
curves and the coupling mechanism. 

We report on a theoretical model utilizing reasonable potential 
energy curves and coupling matrix elements in a two i.tate impact parameter 
calculation. Details are presented in the next section. 



II. THEORY 


We here outline the general theory appropriate to the charge 
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exchange process discussed above. For simplicity v/e discuss only the 
two active electrons. The two electron molecular system for fixed nuclei 
is described by the electronic Hamiltonian 





V +V 


rt 




( 3 ) 


with Ta being the kinetic energy operator for the electrons and the 
potentials V describe the two-body Coulomb interactions amongst the nuclei 
and electrons (See Fig. 1), We treat the colliding system in a quasi- 
classical fashion where the nuclear separation R is determined (in 
principle) from a classical scattering calculation. Therefore H acquires 
an implicit time dependence via R and we need solve a time-dependent 
Schrodinger equation for the wave function of the colliding system. 


14 a •'5' (S'*) 


< ft 3_ ^ (A,t) 

at 



The wave function 
orthonormal basis 



is represented as an expansion of an arbitrary 
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where 


*^ 4 . /#/<&> and b^(t) are the expansion 
coefficients which are to be determined. Various possible choices of the 
bases will be described below. Insertion of Eq. 5 into Eq. 4 yields the 
set of coupled equations. 



In a process where the transition probability is small it may be 
acceptable to truncate the sum to include only the two states of interest 
(we discuss an added complication below). We make this approximation 
calling the incident state i and the final state f. The tv/o basis states 
retained, ^and must have the property that as R -*■ «, they represent 
the separated atom + ion 

However this physical restriction still allows a manifold of choices for 

the bases <f>. Two choices are as follows 

i) Adiabatic bases. 

(a) 

The functions $ are chosen to diagonalize the electronic Hamiltonian 
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(7) 


In this expansion the first coupling term in Eq 6 is seen to vanish and 
transitions are caused by the non-adiabatic coupling term 
We mention that in such a basis set the eigenvalues H^^CR) and (R) 
are unequal for all R unless the states have different symmetry about the 
molecular axis. Therefore H n ^ a ^ (R) can equal hJ^Cr) at the crossing 


mm 


point R v for so-called rotational coupling but not for radial coupling. 


(For radial coupling l I atl S'" while for 

rotational coupling < ^ j wl ’ ere 


©describes the orientation of the molecular axis in a lab-fixed frame.) 


ii) Diabatic Bases. 

Here the functions <f>^ are chosen such that 

Thus the second coupling term vanishes in Eq 6 and transitions are caused 
by the adiabatic interaction 


( 3 ) 
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(d) 


For this basis set H 1 ‘ can equal ' even if n and m have the same 
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symmetry. In general, if the complete expansion, Eq 5, is retained either 
basis set will yield the same results. However, when the basis set is 
truncated one can expect different results depending on the choice made. 

We argue that in a problem of such complexity as this, one can only 
hope to obtain rough estimates of the cross sections and therefore we will 


choose an approximate diabatic basis in which the states <j> are represented 
by the non-interacting atomic wave functions for all R. 


4>c = /* * 

, (10) 
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These bases are approximate since they do not maintain the ortho- 
gonality requirement for all R. Further if the two molecular states have 
different symmetry about the molecular axis one is obliged to use an 
adiabatic representation since H jirn = 0 and transitions are caused by the 
rotational coupling term j 

However rotational coupling is not as important at thermal velocities due 
to its velocity dependence. 

The calculation of the cross section proceeds as follows 
The coupled Eqs6 reduce to ^ 

it <i b - = eK P J /> + 

It 
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with the initial conditions I be (t--oo)j ~ it -~co) s= Q * 
These equations can be solved if the matrix elements H mn are known for 
all R and if R(t) has been obtained. However, in many cases one can 
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obtain an analytic solution to these equations. In particular, if the 
coupling is small then perturbation theory can be used. That the 
coupling is indeed small for Duffendack reactions can be seen by noting 
that the process involves the simultaneous exchange of one electron plus 
the excitation of the remaining electron in a low energy collision. To 
illustrate, in collisions of metal atoms with ionized helium we have 
M = M(ns 2 ), R + = He + lls),M + = M + *(n'acr) t R = He(ls 2 ) where, in the 
incident molecular complex R M the atomic electrons have a zero component 
of angular momentum along the molecular symmetry axis and is therefore 
a l molecular state. Since the coupling vanishes if the projection 
of the orbital angular momentum for the final state was not j , then the 
excited electron must also have a zero projection on the molecular axis. 
Therefore 


he l^ 5 V )/n % l %2) 


where i s given in Eq 3. Since ^ and differ by two orbitals 
the matrix element vanishes identically and must be replaced by a 
second-order coupling determined from perturbation theory 



V„= 2 1 
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the sum on m runs overall virtual, intermediate states. These states 
<{> m are exactly those omitted when the coupled equations 6 were truncated 
to the two state equations 11. The sum over m should be dominated by a 
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single term which represents the configuration in which the ion R + has 
excited the atom into the state M = M (ns, n’to). Thus we take 

4> l = 

<3 *> = ( 

<t> « /i5o;o) /M / / /(r + t ?)> 


The matrix element H . is a multipole excitation of the ns^ orbital 
into the n 1 ^ state by the charge R and as such 
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The constant C . is the multipole moment for the transition in question 
mi 
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The matrix element *^fm describes the exchange of the remaining .ns orbital 
into the 1S R state. For large R it is approximately equal to 

V (|;| 


where I is the energy needed to ionize the ns orbital in M and the 
m 
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integral is a standard overlap integral . From Ref. 8 y/e obtain 
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where 1^ is the ionization energy of hydrogen and a 0 is the Bohr radius, 

(We set a 0 - In = 1 henceforth for notational simplicity.) This result 

is valid when l r « 1„ which is the case here. When this is not the 
f He 9 

case the approximation of Olson may be more appropriate. Finally we 
observe that at large distances the denominator in Eq 13 can be set 
equal to the excitation energy of M 
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Therefore 
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fa ' R 


and the coupled equations for the amplitudes become 

L t ik : = e*p [ “-jr | 


i k = Vft e<pL ~-jr v ’ b, 


The probability for a transition is obtained by setting 

= | b + c+«oj 



where the calculation for(? must be performed over a range of impact 
parameters V yielding the cross section 
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( 22 ) 


Under certain circumstances the coupled equations IT can be integrated 
analytically. Consider the argument of the exponential, which we 


rev/rite as tfct) 

( AR 

4 4*At 


(23) 


At low velocities the denominator of this expression is quite small and 
the integral will rapidly change as the system proceeds along its 
trajectory. The exponential function will therefore oscillate rapidly 
and effectively give no net coupling. However, if the trajectory takes 
the system through a radius R x where Hff(R x ) = Hjj(R x ) [Fig 2] there will 
be a small interval where the integral changes slowly allowing for the 
coupling to be effective. This is the basis of the standard Landau- 
Zener method 6 . However, for very low collision energies an additional 
complication arises. The classical turning point, R 0 , is very close to 
the crossing point R x (see Fig. 2) therefore, the radial velocity R 

cannot be taken to have a constant value as is done in L-Z theory. In 

• • • 

fact R will be changing signs in the vicinity of R x . This complication 
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has been taken into account by several authors however . Since the 


potential curves H^ and H.. 


will cross at large values of R 

A 


for 


Duffendack reactions and since is exponentially small then, one can 
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use perturbation theory 0 to solve the coupled equations 11' . The 

r* 

results^"' are 



Ai is the Airy function as defined in Ref 10 and 


IL c /W* (W* __ 

the reduced mass of the system 


e = = 

energy of the colliding pair 

AF =- 

1 A (llu-ttcc) 1 
i R 1 


fr~ — 

effective raaial force for the 
classical trajectory 


effective potential energy for the 
classical trajectory 

Ha + e£ 

Ri 

and all quantities are evaluated at R = R x . The fact that these results 
are consistent with a v/ave equation analysis if we choose 

(26) 


is discussed in Refs 6 and 7. 
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To proceed further we are required to obtain expressions for and 
which asymptotically take on the values e- z 0 and (See Fig 2). 

Following Nikitin we take for large R 


A,» l «"«-•** 


where 


5 v'zT 
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and 1^ is equal to the ionization energy for the ground state atom M 
(again we have scaled the ionization energies to 1^ and the radius to a 0 ). 
In a similar fashion we have 




with I f equal to the ionization energy of the ion M (n'fl,). That 
and should be exponentially decreasing for large R results from the 
exclusion principle and the observation that the nsa orbitals for the 
incident metal atom are filled as are the Is a orbitals for the outgoing 
rare gas atom. That should decay more rapidly than results from 

..... i 4- 

the more diffuse nature of the M (n's,cr) orbital. In both cases the 
attractive polarization contributions should be negligible at the large 
R considered. The crossing radius R x is then determined by setting 
H... (R x ) = H f f(R x ). Since a.. » we have 


K Cc 
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As discussed above, all of the parameters of the model problem can be 
calculated readily excepting the constants A., A f [or equivalently R 

1 ' A 

and H x = (R x )]. The functions AF, F and V are then given by 


- o^A f t 

t * 1 x +2E ST 


(25') 


V = 

Ri 

where we have assumed aR x » m. Inserting these expressions Into the 
transition probability we can determine its dependence on the impact 
parameter, P , from Eq 24. The P dependence occurs in two places 
U(p) and F(C). Typically the f dependence in F is dropped which permits 
a considerable simplification. However vie retain it for consistency. 
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III. PRELIMINARY RESULTS 


Numerical calculations will be performed shortly. However we can 
qualitatively anticipate the results in two extreme cases. If the argument 
of the Airy function in Eq 24 has a large magnitude we can replace the Airy 
function by its asymptotic value. 
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The transition probability#, can be written as 
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To investigate the circumstances under which the asymptotic forms can be 
used we refer to Eq5 25 1 which yield 


E -V - Etl-ffc) ~k 
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Case A) E » H 

x 

From Fig 2 we see that this corresponds to a small energy defect and 

corresponding large crossing radius R^. The classical distance of closest 

approach R 0 can be less than E and the trajectory will pass through R if 

* x 

p < R x . Under such circumstances we can neglect the contributions to 

2 

from p > R x< If we further replace the rapidly oscillating sin term 
by its average value of 1/2 we obtain from Eqs 22 and 33a 

5rr /.'r t r P, 
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This is the same result as obtained from perturbation theory when the 

6 

proximity of the crossing point to the turning point is neglected. 
Substituting the expressions for and a^ from Eqs 20 and 29 we see that 

as Af 0 (R •* «)^ the cross section is dominated by the exponential terms 
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where I is the energy needed to ionize the ns orbital in M* (ns n'i,) and 
m 

I. is the energy needed to ionize the n'& orbital in Therefore 

T 

I > I- and the cross section decreases as + 0. 
m f T 

Case B) E << H 

From Fig 2 we see that under such circumstances the crossing point can never 
be reached. Therefore can be approximated by Eq 33b for all P . A 
qualitative interpretation can be made. Replacing the phase interference term 
by the average value of 1/2 in Eq 33a we formally have 



Therefore the transition probability when the crossing point cannot be 

(B) 

reached classically, (P^ , equals the transition probability when the 

crossing point is reached, multiplied by a quasi-classical transmission 

coefficient^ describing the probability that the system can tunnel through 
the potential barrier and reach R . An order of magnitude estimate of the 
cross section is then 
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We observe then that as A increases, the crossing radius decreases and the 

coupling at the crossing point increases. Hov/ever since the crossing 

point cannot be reached classically the cross section is dominated by the 

transmission coefficient which decreases rapidly as H grows. Thus the 

X 

cross section for a Duffendack reaction is seen to decrease for a^ *> 0 and 
for Af •* °°. In the former case the decrease is due to a decrease in the 
coupling V.^. while in the latter limit the decrease is due to the inability 
of the classical trajectory to reach the crossing point where transitions 
are most likely. Such a behavior was seen in the experimental work of 
Turner-Smith et. alj 

Two remarks are made in conclusion. If the distance of closest, approach 
R 0 , is much larger than the crossing radius R the transition probability 

A 

5 

given in Eq 24 is no longer correct. In such circumstances Delos and Thorson 
have shown that the correct expression is given by 


<?„ = r^b % £ a-*; 


(3a) 


The quantities x and b are given there and reduce to Eq 24 if the variations of 
V. F and aF are negligible in the interval R 0 to R . Finally we point out 
that to more directly compare the theoretical and experimental work a 
thermal average of the cross sections must be performed, 
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